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Abstract 

Recent results concerning the use of the Correlated Basis Function to investigate the ground 
state properties of medium-heavy doubly magic nuclei with microscopic interactions are presented. 
The calculations have been done by considering a Short-Range Correlation between nucleons. The 
possibility of identifying effects produced by Short-Range Correlations in electromagnetically induced 
phenomena is discussed. 

1 Introduction 

In many-body physics the word correlation indicates effects beyond Mean-Field (MF) theories. In 
nuclear physics it is common to distinguish between short- and long-range correlations. Nuclear collective 
phenomena such as vibrations and rotations are ruled by long-range correlations. These effects arc well 
known and studied since the infancy of nuclear physics. On the contrary the study of the Short-Range 
Correlations (SRC) is a relatively new issue in nuclear physics. These correlations are produced by the 
strong repulsive core of the microscopic nucleon-nucleon interaction at short internucleon distances. In 
spite of the fact that all the microscopic nuclear theories need the SRC, clear signatures of their presence 
in medium-heavy nuclei have not yet been identified. The search for nuclear phenomena showing SRC 
effects is one of the most discussed topics of the last 15 years in the nuclear structure community. 

In this contribution I shall be concerned about SRC. In the first part I shall present some results of 
microscopic calculations regarding ground state properties of medium-heavy nuclei. These calculations 
have been done with realistic nucleon-nucleon interactions and the SRC have been used to tame the 
repulsive core of these interactions. In the second part of the contribution I shall discuss the results of 
a study aimed to identify the presence of SRC in nuclei by using experiments done with electromagnetic 
probes. 



2 The nuclear ground state 



The results presented in this section, have been obtained within the framework of the Correlated Basis 
Function (CBF) theory [Q. The starting point of this theory is the Ritz's variational principle: 
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with the following ansatz for the many-body nuclear state: 
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where |$ > indicates a ground state Slater determinant constructed with a set of orthonormal single 
particle wavefuctions \<f> >, and F is a many-body correlation function defined as: 

i<j \p=l / 

In the above equation S indicates a symmetrizer operator. The two-body correlation functions f p (rij) 
have an operatorial dependence of the same type of that of the nucleon-nucleon interaction. In the present 
context I shall consider correlations up to p — 8 of the type: 

O p=1 ' 8 {i,j) = [1, <r(i) ■ S(i,j)MU) ■ S(<, j)] ® [T(i) • t(j)] (4) 

where S(i,j) is the usual tensor operator. 

The problem is now well defined, and the minimization of the energy functional (Q) is, at least in 
principle, only a technical matter. Monte Carlo techniques have been used || to solve the multidimen- 
sional integrals of cq. ([!]). This approach is exact, i.e. no further approximations are required, but it can 
be used only for relatively light nuclei. So far it has been applied to nuclei up to A=8, and for the 16 
nucleus the Monte Carlo calculations have been done with an approximation consisting in considering 
only up to four-body nucleon clusters ||. 

It is clear that the minimization of the energy functional (|l]) for heavier nuclei requires the use 
of techniques alternatives to the Monte Carlo. These new techniques should simplify the problem by 
including new approximations. The idea is to find approximations affecting only irrelevant parts of the 
calculation, in order to obtain accurate results. 

A technique rather well studied is the Fermi Hypernetted Chain (FHNC) theory [Q, developed first 
for infinite systems such as liquid helium and nuclear matter. A first step in the application of the FHNC 
consists in a topological analysis of the various nucleon clusters formed by the interaction, by the SRC 
and by the antisymmetrization of the many-body wave function. It is possible to sum in a closed form 
all the diagrams with certain topological characteristics; these are called nodal diagrams. The other 
diagrams, named elementary, have to be calculated one by one, as in traditional perturbation theory. 
In this approach the most common approximation consists in neglecting the elementary diagrams. This 
approximation, called FHNC/0, is rather good in nuclear matter, even though it is not accurate enough 
for helium systems since they have a relatively higher density. 

The FHNC theory has been extended to finite nuclear systems in refs. 0, and applied to light doubly 
magic nuclei in refs. [js], ^] for simple interactions and correlations. The main goal of that investigation 
was a test of the validity of the theory. We found that the various sum rules were exhausted at the 0.1% 
level when a specific elementary diagram was included. We named this approximation FHNC-1. 



i« 


4lJ Ca 


full corr exp 
-5.48 -5.41 -7.97 


full corr exp 
-6.97 -6.64 -8.55 



Table 1: Binding energies per nucleon in MeV. The label full indicates that the minimization has been 
done on both correlation and single particle wave functions, while corr indicates that only the correlation 
has been varied. 

The next step of the work consisted in extending the theory in order to use realistic microscopic 
interactions with their full operatorial dependence. This requires to activate the analogous operatorial 
dependence of the correlation function, as given in eq. (|3|). The various terms of the correlation do not 
commute any more with the hamiltonian, and among themselves. The various cluster terms cannot be 
resummed in a closed form 
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Figure 1: Elastic electron scattering cross sections. The full lines have been obtained with the empirical mean 
field wave functions and the V8' interaction. The dashed lines have been obtained by using the fully minimized 
densities with the V8' interaction and the dotted line with the Argonne V14 interaction. 



For this reason it is necessary to introduce in the FHNC-1 computational scheme another approxima- 
tion called Single Operator Chain. The FHNC-1 calculations in this approximation sum, in addition to 
all the nodal diagrams containing the scalar term of the correlation, only those diagrams where a single 
state dependent operator terms appears. This further approximation is slightly worsening the accuracy 
of the result. The sum rules are now exhausted at about the 5% level [0. 

With this approach we have evaluated the ground state properties of the 16 O and 40 Ca nuclei using 
the realistic 2-body interaction Argonne V8' plus the 3-body Urbana IV force J|, |). The first one is a 
reduced version of the Argonne V18 interaction limited to the first 8 channels of eq. ([|). In order to 
compensate for the truncation of 6 operatorial channel the parameters of the three-body Urbana IV force 
have been slightly readjusted 0. 

The minimization procedure implied by eq. ([j]) involves both correlation and single particle wave 
functions. The energies obtained by this complete minimization are compared in tab. ^ with the empirical 
binding energies. The discrepancy may seem rather large, but it is of the same magnitude of that of the 
best microscopic nuclear matter calculations. 

The calculated charge densities distributions of 16 O and 40 Ca have been used to evaluate the elastic 
electron scattering cross section [|| . These cross sections are shown by the dashed lines of fig. [I] where 
they are compared with the data and the results of an analogous calculation done with the Argonne V14 
interaction (dotted lines). Also in this case the agreement with the data is rather poor. Theory and 
experiments start to differ already at 30°. 

The same computational scheme has been used by keeping fixed the set of single particle wave func- 
tions. We used a set already utilized in ref. Q and chosen to reproduce charge radii and single particle 
energies around the Fermi surface. The minimum of the energy functional (Q) is obtained by changing 
only the correlation. The results of this calculation are shown in tab. |l| and in fig. |l| by the full lines. 
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Figure 2: Momentum distribution. The full lines have been calculated considering central and tensor terms of 
the SRC. The dashed dotted lines have been calculated with only the central SRC and the dashed lines show the 
pure IPM results. The white dots are the Monte Carlo results of ref pi. 



The difference between the binding energies obtained with a full minimization and with the one 
restricted to the correlation function only, is of the order of the few percent. This seems to indicate the 
scarce sensitivity of this quantity to the single particle wave functions. On the other hand the charge 
densities, and consequently the cross sections, have been strongly modified. Now the agreement with the 
data extends up to about 60°. 

This last set of single particle wave functions has been used to calculate the momentum distribution 
of the two nuclei considered ||. The result of this calculation is shown in fig. ||. The dashed lines show 
the MF distributions, the dashed dotted lines have been obtained with purely scalar correlations and the 
full lines show the results of the complete calculation. These results show that the contribution of the 
state dependent terms of the correlation is not negligible. The agreement with the Cluster Monte Carlo 
momentum distribution of 16 (white squares) is rather good. 

The large increase of the correlated momentum distributions with respect to the MF results is a 
well known feature produced by the SRC JlCj . Unfortunately the nucleon momentum distribution is not 
directly observable. It is therefore necessary to find observables related to it in order to identify the 
presence of SRC. This search should be done by studying the nuclear excited states. The results of this 
investigation are the subject of the next section. 

3 Search for SRC effects in excitation processes 

In the framework of the CBF theory, the description of the excited states of the nucleus can be done by 
extending the ansatz of eq. (Q): 

|*„>=F|$„> (5) 
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Figure 3: Mayer-like diagrams showing the two- and tree-body clusters terms included in the present work. The 
squares indicate the coordinate where the one-body electromagnetic operator acting on. 



where now |$ n > indicates a Slater determinant containing particle- hole excitations. In principle, in the 
description of the excited states also the correlation function should change, and the variational principle 
(|l|) should be again applied. 

In reality, we have used the approach of Fantoni and Pandharipande which considers both corre- 
lation function and MF potential fixed by the ground state minimization. The transition amplitude from 
the ground to an excited state, induced by an external operator, can be evaluated using cluster expansion 
techniques analogous to those used for the ground state. In the special case of infinite nuclear systems, 
it has been shown ]Tl], [l^] that it is possible to sum all the diagrams of a certain type. 

Our calculations of the electromagnetic responses of finite nuclei, have an additional approximation, 
consisting in retaining only the diagrams that contain a single correlation line. This amount to consider 
only certain classes of two-body and three-body cluster diagrams. As example I show in fig. ^ the 
Mayer-like diagrams considered in the case of an excited state which asymptotically can be described 
in terms of one-particle one-hole excitation. In ref. |l3|] this approximation has been tested against the 
complete calculation for the nuclear matter charge responses. The agreement between the two approaches 
is excellent. The differences are of the order of 10~ 5 . The approach briefly outlined is reliable only to 
describe situations where nuclear collective effects are negligible. This means that the dynamics of the 
excitation should be dominated by single particle excitations. 

For finite nuclear systems, the theory has been applied only for scalar correlation functions. Specifically 
we have used the single particle wave functions and correlations labelled S3 and V8 in ref. || and shown 
in fig. |. 

We first applied our model to high angular momentum stretched states in the discrete spectrum of 
16 O, 48 Ca and 208 Pb |Q. Because of the high angular momentum these states are dominated by a single 
particle excitation. A typical result of this study is shown in fig. [|. The dashed lines represent the MF 
results while the full lines have been obtained by including the correlation terms of fig. 0. Clearly the 
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Figure 4: Two-body correlation functions used in the calculations of the electromagnetic response functions. 

SRC are unable to explain the quenching required to reproduce the experimental data. 

We have applied our model in the quasi-elastic region for inclusive reaction . Also in this case the 
excitation is dominated by one-particle one-hole excitation dynamics, therefore collective phenomena are 
negligible. However, contrary to the previous case, one has to consider the re-scattering of the emitted 
nucleon with the rest nucleus. This effect, usually named Final State Interaction (FSI), reduces the 
maximum values of the responses by a 20-30% factor . 

In fig. U we show the inclusive electromagnetic responses calculated in 40 Ca, without the inclusion 
of the FSI. The inclusion of only the 2-point diagrams of fig. || slightly increases the MF responses (full 
lines in fig. ^J). The additional inclusion of the 3-point diagrams has an opposite effect with respect to 
that of the 2-point diagrams. The final difference between the MF and correlated responses is very small. 
The effect of the FSI is much larger than that of the SRC U]. 

We have extended our model also to describe (e,e'p) coincidence experiments on 16 target JTtJ. To 
show the relevance of the 3-body cluster terms, I show in fig. ^ the normalized differences between the 
one-body response and the calculation considering two- and three body- SRC terms: 



The results of the two calculations presented in fig. [7| have been obtained with the correlations 
functions of fig. |]. Since the two correlations are rather similar, it is plausible to expect that the results 
should not differ very much. This is not the case when only 2-point diagrams are used, as the dashed 
lines of fig. |t] indicate. The inclusion of the three-body clusters, which ensure the correct normalization 
of the nuclear final state, produces results which are rather similar for the two correlation functions. 

The relevance of FSI is clear from the comparison with the 16 O NIKHEF data shown in fig. & 
The full lines have been calculated without FSI and they are distant from the data points, while the 
dashed lines considering the FSI via optical potential are much closer. The SRC are unable to explain 
the need for a spectroscopic factor in order to reproduce the data, as the two lower panels clearly show. 

Our model has been also applied to study nucleon emission induced by real photons. In Fig. ^1 show 
some results of our investigation. The thin full lines have been obtained making a pure MF calculation 
with one-body current only. The inclusion of SRC, specifically the S3 of fig. [|, produces the dotted lines, 
while the inclusion of Meson Exchange Currents (MEC) gives the dashed lines. The thick full lines show 
the sum of all the effects considered. 

The effect of the SRC seems to be relevant at high excitation energies and for large emission angles. 
Unfortunately also the MEC are important in this region, and they produce effects as large as those of 
the SRC, if not even larger. 

We have applied our model to describe (e,e'2p) processes [^(j. In this case the MF responses induced 
by a one-body operator, which in the previous cases were giving the largest contribution to the cross 
section, are not present. 
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Figure 5: Inelastic scattering form factors calculated with SRC, full lines, and pure MF results, dashed lines. 
The left and right panels presents the same results in linear and logaritmic scale. 



In fig. |l0| the results for the 16 (e,e'2p) under the kinematics explained in the caption are shown. 
The dotted lines represent the results obtained by using only two-point diagrams, the dashed ones by 
adding the tree-body diagrams and the full lines include also the contributions of the excitation of the 
A resonance. It is interesting to notice the changes of the relative importance of the various components 
of the calculation by changing the nuclear final state. The common features of all the results shown in 
the figure is that the 3-point diagrams reduce the effect of the 2-point ones. This effect of the 3-point 
diagrams is consistent with all our previous calculations, and its physical origin is related to the correct 
normalization of the nuclear final state, as it has been detailed discussed in ref. |]l5| . 

Our calculation show that, even in this case, the contribution of the two-body currents, consisting 
here only in the A excitation term, may cover the SRC effects. The A excitation dominates the (lpl/2)~ 2 
+ and the (lpl/2) _1 (Ip3/2) _1 1 + cross sections. In the other cases considered its effect is relatively 
small. 



4 Conculsions 

Microscopic calculations of nuclear properties require SRC to tame the strong repulsive core of the 
realistic nucleon-nucleon interaction. In our theory the SRC are an input fixed by the minimization of 
the hamiltonian expectation value. Within our theory we succeeded in calculating the properties of 16 
and also of 40 Ca. The discrepancy with the empirical values of the binding energy is of the same amount 
of that of the best nuclear matter calculations. The value of this quantity shows a scarce sensitivity to 
the change of the single particle wave functions, contrary to the density. 

To search for effects clearly produced by the SRC we developed a model describing nuclear responses. 
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Figure 6: Charge, Rl, and current, Rt, inclusive responses. Full lines MF responses, dotted lines MF plus 
2-point diagrams, dashed lines MF plus all the diagrams of fig. ^. 
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Figure 7: Normalized differences, eq. ^| calculated with the correlation functions of fig. ||. The dashed lines have 
been calculated by using two-point diagrams only, while the full lines show the inclusion of tree-point diagrams. 




Figure 8: ls O (e,e'p) 15 N reduced cross section for the emission of the lpl/2 1 proton calculated with the 
correlation functions of fig. ^. The full lines show the results obtained with a purely real MF basis, the dashed 
lines have been obtained with an optical potential. The curves in the lower panels have been produced by 
multiplying the upper curves with spectroscopic factor of 0.7. 
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Figure 9: 16 O (7,p) cross sections as a function of the nucleon emission angle for various values of the photon 
energies. The thin full lines show the MF results. The dotted lines include the effects of the S3 correlation. The 
dashed lines the MEC and the thick full lines all the effects. 
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Figure 10: 16 (e,e'2p) 16 cross sections as a function of the emission angle of one proton. The initial energy 
of the electron is 500 MeV, the nuclear excitation energy 200 MeV, the momentum transfer 300 MeV/c, and the 
second proton is emitted at 60 . The various labels in the figures indicates the single particle levels where the 
two protons are emitted and the angular momentum and parity of the residual A-2 nucleus. The dotted lines 
have been obtained with 2-point diagrams only, the dashed lines by adding the 3-point diagrams and the full lines 
consider also the contribution of the A excitation. 
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Within this model we have investigated the role of the SRC in (e,e') experiments both in the discrete and 
in the continuum spectrum. We have also investigated the single nucleon emission with both real and 
virtual photons. These processes are dominated by the MF transition densities, and the SRC produce 
very small effects. In specific situations when they show up, also the two-body MEC contribute, and the 
effects of these two processes strongly mix. 

The two-nucleon emission processes are more promising, because in this case the MF contribution 
is absent. Also in this case the SRC effects fight against the MEC ones produced by the A excitation. 
On the other hand, the large number of variables could allow for the possibility of disentangling the two 
processes. 
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